Entropy-induced Microphase Separation in Hard Diblock Copolymers 
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Whereas entropy can induce phase behavior that is as rich as seen in energetic systems, mi- 
crophase separation remains a very rare phenomenon in entropic systems. In this paper, we present 
a density functional approach to study the possibility of entropy-driven microphase separation in 
diblock copolymers. Our model system consists of copolymers composed of freely-jointed slender 
hard rods. The two types of monomeric segments have comparable lengths, but a significantly 
different diameter, the latter difference providing the driving force for the phase separation. At 
the same time these systems can also exhibit liquid crystalline phases. We treat this system in 
the appropriate generalization of the Onsager approximation to chain-like particles. Using a linear 
stability (bifurcation) analysis, we analytically determine the onset of the microseparated and the 
nematic phases for long chains. We find that for very long chains the microseparated phase always 
preempts the nematic. In the limit of infinitely long chains, the correlations within the chain become 
Gaussian and the approach becomes exact. This allows us to define a Gaussian limit in which the 
theory strongly simplifies and the competition between microphase separation and liquid crystal 
formation can be studied essentially analytically. Our main results are phase diagrams as a function 
of the effective diameter difference, the segment composition and the length ratio of the segments. 
We also determine the amplitude of the positional order as a function of position along the chain at 
the onset of the microphase separation instability. Finally, we give suggestions as to how this type 
of entropy-induced microphase separation could be observed experimentally. 



I. INTRODUCTION 



Microphase separation (MPS) is the phenomenon 
where an initially homogeneous phase develops an inho- 
mogeneous spatial structure on a microscopic scale. Usu- 
ally such systems consist in part of thermodynamically 
incompatible components that left by themselves would 
tend to (macroscopically) phase separate. However, due 
to additional constraints of a physical or chemical na- 
ture the spatial separation between the components is 
prevented from increasing beyond a microscopic length 
scale. This leads to phases in which the components can 
demix only locally. There are a few archetypical examples 
of systems showing MPS: (i) Two (usually flexible) poly- 
mers species that have an unfavourable mutual interac- 
tion energy which are joined together by a chemical bond. 
This type of block copolymers [E S 01 shows a wealth of 
microphases. (ii) Side-chain liquid crystalline polymers 
(LCPs). Here, liquid crystal- forming groups are linked to 
polymer backbones through flexible spacers. The most 
prominent phase of these systems is the smectic, where 
the LC groups form orientationally ordered layers sep- 
arated by disordered lamellae containing the poymeric 
backbones 0, Q . (iii) Ternary systems consisting of wa- 
ter, oil and an amphiphilic surfactant. These systems 
can show a variety of microstructured phases, with the 
amphiphilic surfactant stabilizing the oil-water interfaces 
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and thus preventing "macrophase separation" Q . All 
three of the cases above are examples of thermotropic 
systems, i.e. systems in which the phase behavior is gov- 
erned by temperature as a controlling variable, reflecting 
the dominance of energetic effects. 

Recently, MPS was observed in an entirely new class of 
systems. Binary mixtures of bacteriophage viruses and 
(small) latex spheres with varying size ratios showed a 
surprisingly rich phase behaviour, including a lamellar 
phase 8J. In this phase, the lamellae are defined by a 
"smectic" arrangement of the rodlike virus particles in 
layers with the spherical latex particles in between the 
layers [113. These results are remarkable for two rea- 
sons. First, unlike the previous archetypal cases of MPS, 
we are dealing with a binary mixture which phase sep- 
arates on a microscopic scale. There is no "hard" con- 
straint like a chemical bond that prevents the two species 
from phase separating on a macroscopic scale, and both 
species remain in a fluid state within the layers. Second, 
it was argued that the virus particles as well as the la- 
tex spheres can be modeled to a good extent to interact 
as hard bodies. Consequently, the driving force causing 
this MPS must be of an entropic nature. This is also 
in stark contrast with MPS in block copolymers, LCP's 
and amphiphiles where the dependence on temperature 
is strong and hence indicates a predominantly energetic 
effect. The possibility of this type of MPS was already 
explored in computer simulations 10] and found to be 
qualitatively well described within the so-called second 
virial approximation [Tol lllj | , the validity of which can 
only be guaranteed at low densities. However, as the ex- 
perimental systems are far from dilute this last treatment 
may not capture all the essential ingredients. It has been 
argued that MPS in binary mixtures may be caused by 
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the so-called depletion effect H, , which is generically 
a many-body interaction effect and is not well described 
with a second-virial theory. Consequently, a more ac- 
curate approach would be required, certainly in order to 
resolve in detail what prevents the system from demixing 
macroscopically. 

That entropy per se can be the driving force for phase 
transitions has by now been well established. There 
are many examples ranging from ordering in monodis- 
perse systems like the liquid-to-crystal transition in hard 
spheres [T^ and the isotropic-to-nematic transition in 
slender hard rods lll^ . to demixing in binary mixtures, 
like e.g. the Asakura-Oosawa (AO) mixture of hard 
spheres and ideal spheres, which is used as a model for 
colloid-polymer mixtures [l^ ITEl IT6| . In essence, the 
physical mechanism in all these systems is the same; the 
gain in effective "free volume" available to the particles 
upon ordering offsets the loss of entropy of disorder or 
mixing respectively. For the AO mixture this is usually 
referred to as the previously mentioned depletion effect; 
the ideal polymers are depleted from a shell around the 
impenetrable colloids. Overlap of these depletion shells 
increases the free volume available to the polymers and 
hence this system phase separates into a colloid-rich and 
a colloid-poor fluid 17J. However, whereas entropy can 
induce phase behaviour that is at least as rich as seen in 
energetic systems, MPS remains a very rare phenomenon 
in entropic systems 18]. 

A variant of the depiction effect was also recently dis- 
covered in theoretical treatments of binary mixtures of 
thin and thick hard rods 19, 2^- These systems are 
seen to be able to phase separate in two isotropic fluid 
phases due to depletion. Here, however, the depletion 
interaction appears as a genuine two-body effect 20], in 
contrast to the AO system, in which it is predominantly 
a three-body effect. Consequently, this form of depletion 
effect survives the Onsager limit (length^^width) applied 
to both species of rods, and for sufficiently asymmetric 
rods, preempts the usual transition to the orientationally 
ordered nematic phase |3 US 1^ • These predictions 
have since been corroborated by simulations . In 

the present paper, we propose to use the two-body de- 
pletion effects between slender rods of different diameters 
to construct a system which shows entropy-induced MPS. 
Taking our cue from the concepts developed in the field 
of thermotropic block copolymers, we connect a chain 
of freely rotating "thick" hard rods to a chain of freely 
rotating "thin" hard rods. The above-mentioned un- 
favourable depletion interaction between these two types 
of rods provides the tendency to fully demix, whereas the 
joint (connecting the two strands) prevents this. The so- 
constructed system of freely jointed hard diblock copoly- 
mers (HDC) is in our view one of the most simple systems 
conceivable showing entropy-induced MPS. Furthermore, 
and contrary to the case of MPS in the binary rod-sphere 
mixtures, the physical mechanism is both clear and ro- 
bust. Of course, there is as yet no direct candidate for 
an experimental system well described. However, it may 



certainly be possible for experimentalists to connect (pos- 
sibly long and flexible) chemically inert polymers to the 
ends of virus particles like TMV [Ull^. Together with 
an appropriate solvent this may mimick an effective rod- 
coil system with only hard body interactions. In this 
system, the polymer tails are likely to stabilize the smec- 
tic phase of the virus particles and this could be viewed 
as a microseparated phase. 

In order to describe this system we employ a den- 
sity functional theory in the second-virial or Onsager ap- 
proximation starting from first principles. We assume 
that multiple overlaps between two chains as well as self- 
overlaps of the chains are unimportant. All three of the 
above ap pro ximations, common in theoretical treatments 
of LCPs ^mmUg, should become exact in the On- 
sager limit where the lengths of the rods involved is much 
larger than their widths. The stationarity equations that 
determine the stable phases in our theory are solved lo- 
cally by means of a bifurcation (or, equivalentl y, l inear 
stability) analysis of the isotropic fluid phase Is^ . 
Apart from fluctuations with a nonzero wave vector cor- 
responding to a microseparated phase, we also consider 
spatially homogeneous fluctuations with nematic symme- 
try, in order to study the competition between these two 
types of ordering. For both phases, we obtain closed ana- 
lytical expressions for the spinodal density. We find that 
for long chains and nonzero difference in the widths, the 
microseparated phase always preempts the nematic. 

Naturally we want to make contact with the vast 
amount of literature on thermotropic block copolymers 
in the weak segregation limit. Most of these follow the 
original treatment proposed in the seminal paper by 
Leibler [s^. Leibler considered diblock copolymers in- 
teracting via the heuristic Flory parameter x con- 
structed a Landau expansion in the average composition 
fluctuations. By applying the "random phase approxi- 
mation" and retaining only leading orders of the Fourier 
modes, he was able to map out more or less the complete 
phase diagram. Subsequent refinements extended the 
theory to the strong segregation regime l34l , added fluc- 
tuations [35| and included extra phases |36| . but did not 
change the essence of the approach. Leibler's results have 
been confirmed qualitatively by experiments (Ref. P| and 
Refs. therein) a nd, for finite chains lengths |3q ]. by sim- 
ulations (Refs. 113,133 and Refs. therein). The Leibler 
approach treats the correlations within the polymers on 
the Gaussian level We can therefore connect to this 
approach by applying the Gaussian limit to our model 
of freely jointed HDC's. Within this limit our theory 
becomes equivalent to that of Leibler as far as the treat- 
ment of the intracham interactions is concerned. How- 
ever, the mterchain interactions between the polymers 
are essentially different in the present case, as they are 
of a geometric nature, i.e. totally fixed by the dimensions 
of the composing hard rods. In the Leibler theory, these 
interactions are described generically by means of the 
freely adjustable Flory parameter. A full exploration of 
the parallels between the two approaches, however, was 
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beyond the scope of this work. 

Another class of systems, that appears as a special 
case of our model are the well-studied rod-coil diblock 
copolymers. These consist of one stiff (rodlike) block 
and a much more flexible part. In such systems, liq- 
uid crystalline ordering competes with MPS and a num- 
ber of theoretical studies have been devoted to the 
subject. Most of these combine the Leibler approach 
with an additional Maier-Saupe anisotropic orientational 
interaction resulting in the appearance of a nematic 
phase (and sometimes an additional smectic phase) in 
the phase diagram, besides the various microseparated 
phases [H Hfl El El El H El . However, the ratio 
of the Flory and the Maier-Saupe interaction parameters 
in these approaches is rather arbitrary, whereas in the 
present approach microseparated and nematic ordering 
both result from the same geometric origin with no room 
for additional tuning. 

Finally, there has been some related work on more ide- 
alized but conceptually simpler systems in the context of 
entropic liquid crystals. Holyst considered parallel nail- 
shaped particles which showed a nematic-to-smectic Ad 
transition 46]. As a model for surfactants, Bolhuis and 
Frenkel studied non-additive complexes of hard spheres 
and ideal spherocylinder-tails 42] where Schmidt and 
von Ferber used hard slender rods for the tails E^ Of 
particular relevance to the present work is Ref. l49l] where 
Diichs and Sullivan investigate the phase behavior of hard 
wormlike diblock copolymers. However, in this latter 
work only differences in persistence length are consid- 
ered and not in thickness between the two components. 
Consequently they only find competition between a ne- 
matic and a (orientationally ordered) smectic phase, in- 
stead of the (orientationally disordered) lamellar phase. 
Moreover, only numerical solutions to the stationarity 
equations are presented, whereas we are able to obtain 
additional analytical insight through the stability analy- 
sis of the isotropic fluid phase. Lastly, van Duijneveldt 
and Allen used Monte Carlo simulations to study the ef- 
fect of flexible tails on the phase behavior of spherocylin- 
ders 111] . This was later extended by Casey and Harrow- 
ell to rod-coil molecules of which the isolated rods do not 
posses a smectic phase [sil ]. 

Although our theory is formulated for chains with a 
finite number of rodlike segments, we devote the major 
part of this paper to chains with an infinite number of 
segments in which the correlations between the segments 
are Gaussian. We formulate a consistent Gaussian limit, 
in which the number of model parameters reduces to just 
three. The limit is chosen in such a way that we can 
still consider the competition between MPS and nematic 
ordering. This is achieved by letting the difference in 
thickness between the two types of rods to become in- 
finitesimally small. The limit moreover is such that most 
of the assumptions in the original derivation of the model 
are fully satisfied. The most prominent results are phase 
diagrams as a function of the model parameters, showing 
the regions of stability of the microseparated or nematic 




FIG. 1: An example of a hard diblock copolymer. A freely- 
jointed chain of Ma hard rods with dimensions I a and cLa (left 
side) are connected to a freely-jointed chain of Mb hard rods 
with dimensions and de (right side). 



phases. Furthermore, exploiting the features of the bifur- 
cation analysis, we are able to calculate the relative order 
along the polymer in the microseparated phase at the bi- 
furcation point. The outline of the paper is as follows: in 
Sec. H we define the model and develop the formalism. 
In Sec. HI we briefly discuss the symmetry of the phases 
involved. The bifurcation analysis is the topic of Sec. IV 
and the Gaussian limit is applied in Sec. V. Sec. VI is the 
results section and we end with a discussion in Sec. VII. 



II. MODEL AND FORMALISM 

We consider a monodispersc fluid of N diblock copoly- 
mers in a volume V . Each polymer is a chain of freely- 
jointed cylindrical rods connected end-to-end where the 
first Ma rods are of type t = A having length /a and 
width c^A and the last il/e rods are of type B with di- 
mensions l-Q and o?B (see Fig. We assume that both 
types of rods are very slender, 1^ ^ dr, with r e {A, B}, 
hard bodies, i.e. impenetrable to other rods. The to- 
tal number of segments in a chain is M = Ma -I- Mb 
and every segment has a label m S {!,..., M} specify- 
ing its position in the chain. The state of a segment is 
described by the position of its center of mass and 
an orientation, being a unit vector tUm pointing along its 
long axis in the direction of increasing m. The config- 
uration of a whole chain ^ is fully characterized by the 
position of one of its segments (say the first; ri) and 
the orientations of all of them, SI — {lui, . . . ,ujm}, so 
^ = {ri, ft}. The position of a segment m is then given 

by = ri -I- i X]r=~i i^kt^k + h+ii^k+i) where h = I a if 
k< Ma and lBiik> Ma + 1. 

In density functional theory (DFT) the free energy of a 
(possibly inhomogeneous) fluid of molecules is expressed 
as a functional of the sin gle- molecule configuration distri- 
bution function, p^^^{^) jH]. Using the second- virialjor 
Onsager) approximation it is formulated as follows |5,1| 



= j d^p^'H^) [log(VTp(^)(0) -1 

The integrals are over single-molecule configuration space 
where J d^ = J drodfl and J dfl = J Yim doJm and 
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J diu = Jq^ dcj) dOsinO. Further, p^^\^) is normalized 
as follows / p'-^H^)d^ = N. The factor (3 equals {keT)-^ 
in which is Boltzniann's constant and T the temper- 
ature. The volume Vt we call the 'thermal volume' and 
is a product of the de Broglie thermal wavelengths of the 
molecules jHllil. The quantity <^{^,^') is the Mayer 
function of two molecules with configurations ^ and 
As we are dealing with hard segments, the potential en- 
ergy between two chains is cxi when they overlap 
and when they don't. Consequently, the Mayer func- 
tion is given by 



— 1 if overlap 
if no overlap 



. (2) 

The configurations of both chains involved can be highly 
irregular and the dependence of $ very complicated. 
Therefore we approximate the chain-chain Mayer func- 
tion $ by the sum of all the segment-segment Mayer 
functions 4>m,m', 



M 

E 



)■ (3) 



This expression, to which only individual pairs of seg- 
ments contribute, is actually the first term in a system- 
atic expansion of the Mayer function. Higher order terms 
involve interactions between more than two segments si- 
multaneously 1.30] . At this point we note that apart from 
(i) the second virial approximation and (ii) the above 
expression for the chain-chain Mayer function, another 
(iii) approximation has been made. In this formalism 
the chains are allowed to self overlap, i.e. other than 
the spatial constraint that successive segments are con- 
nected to each other there are no interactions within 
the chain. All three of these approximations are com- 
monly used and corrections to the first two are small 
when > d,- nil. The neglect of the effects 

of self-overlap is assumed to be reasonable in a dense 
polymer melt [s^ where screening effectively compen- 
sates the intramolecular interactions and as a result in- 
teractions between distant parts of the same chain are 
indistinguishable from interactions with the average en- 
vironment because of loss of intrachain correlations. 

In thermodynamic equilibrium, the free energy reaches 
a minimum and the functional is stationary. Therefore, 
we consider the variation of Eq. Qwith respect to p'^^\ 



(3T -P^x^Q 



(4) 



with the chemical potential p, playing the role of Lagrange 
multiplier needed to enforce normalization. Eliminating 
p from Eq.Qlyields the (selfconsistent) stationarity equa- 
tion. 



iVexp [/d$V^)(O<l>A./(eg0] 
/d|exp[/d^V(iHl')*A/(^,^'): 



(5) 



In order to proceed, we define the single-segment dis- 
tribution function (SDF) (of segment m), pm{Ym,^m), in 
the following way 



p„i(r„i,w,„) ^ j W dukp'^^Hi) 

= J II doJkP^^Hri{rra,n),n), (6) 



in which ri is given by r,„ - i YlT=ii^ktO'<: + h+itOk+i) 
and the product is over all segments k but the mth. Inte- 
grating Eq. ISjover all (life except for (2)„i as well and using 
Eq. |2|we obtain a set of equations, 

N f -i-r 

M „ 

/ d4,dt:'fc,pfe'(r'fc,,(i^,)'/'fe,fc'(rfe -r'fe,,c2'fe,w^,) 

,k'=i'' 



where Q is the normalization factor; i.e. the SDF's are 
normalized in the same way as p^^^ : i.e. / drdujpmiv, uj) — 
N. 



III. PHASE BEHAVIOUR AND ORDER 
PARAMETERS 

A. Isotropic Phase 

At low polymer number density, n = N/V, the sys- 
tem is in the isotropic fluid phase, and pm{vm,i^m) is a 

constant, so due to normalization, pm°^ = n/An. Conse- 
quently, 



dr'dLo'pi'°Uk,k'{vk-r\u;k,to') 



(7) 



dio'lklk' {dk + du' ) sin 7(a)fc • lj') 

(8) 

= ~\Tmlklk'{dk +dk'), 

where 7(0) • tD') is the planar angle between Co and uj' 
and one can recognize lk^k'{dk + dki)sai^{il),Q') as the 
excluded volume of two rods k and k' with respective 
orientations uj and uj' . This yields the following normal- 
ization factor in the isotropic phase 



(47r: 



Ma + Mb 



exp [-iym {Mllj^df^ + Mf^M^lp^l^iidf^ + de) + Mllld^)] 

Choosing the dimensions of rod A as units, we define 

I^Ib/Ia, d^ds/dA, Mb/Ma, (10) 



(9) 
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and a dimensionless segment density in a symmetric way, 



Then, Eq.^ becomes 

{I + Ml){l + Mid) 



(11) 



4 (l + Af)(l + M/2d) 



(12) 

where we have also used M = Ma + M-q ■ We also note 
that the normalization factor Qiso is exactly the partition 
sum of the block copolymers in the isotropic phase. 



B. Nematic Phase 

In the (uniaxial) nematic phase, there is orientational 
order with respect to a direction n, however, the system 
is still spatially homogeneous. Therefore, the SDF can 
be expanded in Legendre polynomials. 



(r , a)) = p„ (cj) = n V — a^l^ Pj (w • n) 



with coefficients 



,0) - 



dCj' Pj{ili' ■ n)/9,„ ((!;'). 



(13) 



(14) 



Due to normalization, a^m = 1 as it is in the isotropic 
phase and due to up-down symmetry of the nematic, all 



for odd J (in the isotropic fluid, a 



U) 



for 



all j 7^ 0). The lowest-order coefficient different in the 

(2) 

nematic and the isotropic phase is am which is the usual 
Maier-Saupe order parameter. The physical incentive 
to form a nematic is that the average excluded volume 
between rods is smaller (and therefore the average free 
volume available to the rods is larger) in the nematic 
phase 0. 



C. Microseparated Phase 

Microseparated phases consist of spatially distributed 
regions rich either in type-A or type-B rods and are typ- 
ically governed by a single dominant wavelength. These 
phases exist in a variety of types exhibiting various de- 
grees of symmetry, e.g. lamellar, hexagonal, bcc and even 
more exotic morphologies like the gyroid phase P, [s^ . 
In this paper we do not consider the various symmetries 
of microseparated phases but focus on the magnitude of 
the dominant wavelength and the density for which it be- 
comes unstable. To that end, we observe that the SDF 
can be expanded in terms of plane waves 



Pm(r,'^) 



qGL 



(q,^)e'' 



(15) 



with L some set of wave vectors and the "coefficients" 
given by 

p„(q,cD) = ( dv'e-'^-^' pm{v',uj). (16) 



In general there will be orientational order within the 
domains and consequently the coefficients still depend 
on the orientation. If needed, one could proceed and 
expand these coefficients again in spherical harmonics. 
However, in order to simplify the analysis, this additional 
orientational order in the microseparated phase is usually 
neglected which we will show in Sec. is permitted for 
the case of infinitely long polymers. In homogeneous fluid 
phases like the nematic, the SDF is independent on the 
spatial coordinate and only the coefficient /5„i(0,a)) at 
zero wavenumber survives. 



IV. BIFURCATION ANALYSIS 
A. The Bifurcation Equation 

At low densities, the isotropic phase is the globally sta- 
ble phase, but at higher densities it will become unstable 
with respect to lower symmetry phases exhibiting some 
form of ordering. Points where these lower-symmetry 
solutions branch off the isotropic solution are called bi- 
furcation points and the densities at which this happens, 
bifurcation densities. Different solutions may bifurcate at 
different densities from the isotropic phase. Generically 
the particular solution which bifurcates at the lowest den- 
sity, will give rise to the first ordered phase that is also 
thermodynamically stable with respect to the isotropic 
phase. In this section, we perform a linear stability (or 
bifurcation) analysis around the iso trop ic parent solu- 
tion, along the lines of Refs. |3(J. ISlL l3^. Consequently, 
we assume isotropic distributions with a perturbation of 
lower symmetry. 



Pm{r,oj) = 1- ep,„,i(r,tj), 

in 



(17) 



where the proper normalisation of the SDF requires 
J drdujpm,i{j^,uj) — 0. Inserting this in the stationar- 
ity equations [71 we linearize the exponent with respect to 
the infinitesimal parameter e. 



exp 



^ / dr'du;'pk'{r',Lu')(l)k,k'{rk~r',uJk,u}') 

k,k' = l '' 

M 

\'Kn ^ lkh'{dk + dk') 



exp 



k,k> = l 



k.k' 



1 + e^ / dr' du;' pk', i{r', Lu')(t)k,k' (r k -'r',LUk, to') 



(18) 
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Equating orders in e, to zeroth order, we re-obtain the 
isotropic result, Eq. O To first order this yields the so- 
called bifurcation equations, 



^ / dr'dw'pfe',i(r',w')0/c,fc'(r/c - (19) 

k,k' •' 

These can be interpreted a generalized linear eigenvalue 
problem with eigenfunctions /9m,i(r, cD) and eigenvalue n, 
the bifurcation density. There is an infinite hierarchy 
of solutions to Eq. 1191 for varying degrees of symmetry. 
However, we are only interested in the one (or the few) 
corresponding to the lowest bifurcation density. Note 
that the explicit dependence on the normalization factor 
Q has dropped out since integration over r,„ and 
trivially yields zero on the left hand side by definition 
and, after rearrangement of the integrals made possible 
by the finite range of the Mayer functions (fik.k' j also on 
the right hand side. 

In order to make the bifurcation equation. Eg. 1191 more 
transparent we define for the moment as an auxiliary 
quantity, the fields, 

(20) 

in terms of which the bifurcation equation becomes 



Pm,l(rm, <JJm) 



^ ' k " k"yim 



J dujk"Hk{rk, 



i^k 



(21) 

However, this field Hk is a function of rfe and ujk whereas 
on the left of Eq. 1211 we have a function of and LOm- 
And these are not independent and as 



rm = rfe + Vk.r: 



(22) 



where the vectorial "path" Vk m between k and m is given 
by 

m — 1 

'Pk,m — 5 ^ {Ik'LOk' + h' + l'^k' + l) , (23) 
k'=k 

for k < m. Further, Pm,m = and the case of fc > m 
can be obtained by realizing that Vk,m = —'Pm.k- Conse- 
quently, the interlying orientional integrations in Eq. 1211 
have to make the connection and "transfer" the field from 
segments k to m. 

We return to Eg. 1191 and insert the constraint. Eg. 1221 
via a delta function 



Pm,l(rm, IjJm) 

n \ ^ 



(4ttW ^ _ 

^ ' k " k"=jtm 



j J\ d,Lbk" J drkS{r„i - Tk +V,n,k) 



Next, we observe that in Eq[^ there appear two spa- 
tial convolution integrals. Therefore, it makes sense 
to proceed with a Fourier transform (i.e. ^(q) = 
V^^ J dr,„e"*i '"'"5(r,„)), yielding 



Pm,l(q, Wm) 



(47r 



k k" 



^ / duj'pk',l{q,Ld')^k,k'{(l,'^k,'^')- (25) 
k' 

This is the general form of the bifurcation equation for 
a lower-symmetry solution bifurcating off the isotropic- 
fluid parent solution. Note that the q- vector is the same 
for all segments. Furthermore, at this point, we have not 
yet specified the internal structure of the polymer, only 
that it is a chain of cylindrically-symmetric (rodlikc) seg- 
ments which contains no closed loops. Concerning the 
rodlike segments, the Fourier transformed Mayer func- 
tion (f)k.k' is calculated in Appendix ^ and is for very 
slender segments [Ik S> dk) given by, 

0fe,fe'(q,(ife,a)fe-') = -Ikh'idk + dk')\tOk x '^fe'|x 

jo{^lkCl-i^k)jo{^lk'<i-<^k'), (26) 

where we have already discarded higher-order terms con- 
taining (dk + dk')c{ as the wave vector will be at most of 
order 1/Ia,b so these terms will be small. The function 
jo{x) = sinx/x is the spherical Bessel function of zeroth 
order. We proceed by solving Eq. |2Hlto which we refer 
as the bifurcation equation from now on. 



B. Nematic Solution 

We first consider the nematic solution, which is also 
the simplest. In the nematic phase, q = and the orien- 
tational integrals in the bifurcation equation are trivial 
and it reduces to 

Pm,l('^m) = J- X! d(2;'pk'.l{<^')(i)7n,k'{t^m,<^'), (27) 
^ fc' 

where PmAi^m) = /5m,i(0,Wm) and 

(t>m,k'{'^m,'^k') — —Imh' {dm + dk')\LUm X a)fe/|. (28) 

is simply minus the excluded volume of two rods with 
fixed orientations, ujm and ujk' ■ This bifurcation equation 
is the same as that of a mixture disconnected rods j2£|, 
so for orientational ordering the connectivity of the rods 
within the chains does not play a role. The kernel (f>m,k' 
is now only a function of the planar angle 7 between 
the orientations of the rods, jwm x Lb'f.,\ = |sin7(a;m • 
(2>k')\- Consequently, due to this uniaxial symmetry the 
eigenfunctions of 4>rn,k' and therefore of Eq.[57|are simply 
the Legendre polynomials Pj (see Appendix IB)| . 



J2 J dr'dLo'pk',i{r',Lu')(bk,k'{rk-r',Lbk,Ld'). (24) / 



dLu'^rn,k'ito,to')Pj{'^' -n) = -Imh' (dm + dk')Sj Pj {io -n) , 

(29) 
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with Sj the Legendre coefBcients of | sin7|. In case of the 
nematic phase, it is well known that this becomes first 
unstable with respect to the mode j — 2, so i(a),„) = 

{5/4:TT)ncm' P2{(jJm ' n) with Cm the Legendre coefficients. 
Then, the bifurcation equation becomes 



„(2) 



k' 



(30) 



with S2 — — 7r^/8. This is an M x M matrix eigen- 
value equation and therefore in principle much too large 
to solve. However, by observing that the geometric fac- 
tor on the right hand side does not so much depend on 
the segments m, k' but on whether they belong to A or 
B, we can split the summation, ^j., = X^rSfe'er with 
T = A, B. Then, we can define the "type-average" co- 



efficients, Cr^"* 



Emer Cm and Eq. EOl bccomcs. 



t'=A,B 



(31) 



Rewriting this in terms of dimcnsionlcss quantities. 



C2 = 



32(1 + MPd) 



(32) 



with 



1 \MI{1 + d) 
1{1 + d) MPd 



and 02 



„(2) 

-A 

„(2) 



(33) 

we now have reduced the problem to a simple 2x2 ma- 
trix eigenvalue equation. There are two solutions for the 
density, 



32(1 + MPd) 
n± = — ^ ^x 

TT 

(trGs ± \/tr2G2 -4det G2) /(2detG2), (34) 

with det and tr denoting the determinant and trace re- 
spectively. As the determinant of G2 is negative, only the 
minus sign in Ea. l34l vields a positive bifurcation density 



_ 32{1 + Ml'^d) 

^ncm — ^ 
IT 

(trG2 - Vtr2G2 -4det G2) /(2detG2). (35) 

Note that, within the context of the model as introduced 
in Sec. ^ this analytic expression for the nematic bifur- 
cation is an exact result. In the wider context of liquid 
crystalline polymers, a more general derivation of the ne- 
matic bifurcation density can be found in Ref. |32j |. 



C. Microseparated Solution 

In a microseparated phase, the wave vector q is 
nonzero and the orientational integrals in the bifurca- 
tion equation have to be performed explicitly. However, 
we can make much progress by observing that most of 
the integrals are still trivial, i.e. if segment k" does not 
lie between k and m it does not help to "pass on" the 
infinitesimal field Hk or equivalently, there is no depen- 
dence in the factor exp(— iq • Vk,m)- Consequently, these 
M — \m — k\ — 1 integrations each contribute a factor 
^ dCo ^ At: which is in total (47r)A^-l™-'=l-i. On the 
other hand, concerning the intermediate segments k" be- 
tween k and m; the only dependence on ujk" is in the path 
'Pk.m- Therefore, suppose for a moment that fc + 1 < m. 




\k"=k 

and it is easy to show that 



(36) 



dcDe-'i-''="" = 47r 



sin qk 

qh" 



^4TTjo{qlk"), (37) 



where we have used q = qq with q being the length and 
the unit vector q the direction of the wave vector. When 
m+l < fc, there is an extra minus sign as Vm,k = —'Pk,m 
but this does not change the result |^ only the end fac- 
tors in Eq. 1201 Consequently, we define the factor 



Fk,'m{q) 



n™'Jfc+i Jo(g'fc") for fc < TO- 1 



1 



for fc = TO — 1 , TO 



(38) 

which is symmetric so Fk^miq) = ^m.fc('z)- Inserting this 
in the bifurcation equation yields, 

/ (it2)'pfc',i(q,a)')0,„,fc/(q,(D„i,(D')+ 



fc' 
n 



(47r) 



k^jn k' 

dujdw' cos [\ci-lkLo) pk'.i{<l,'^')^k,k'{(l,<^,<^'), (39) 



where ak,m = sign(TO — fc) is the sign of to — fc. Instead of 
the other "end factor" exp (um.fc^jq ■ hi^) we have used 
cos (^q • Zfcw) as within the integral only the even part 
in q survives. The first term on the right hand side is 
due to the infinitesimal field Hm directly on segment to; 
the second term contains the contributions Hk on seg- 
ments k ^ m which are being transferred to segment to 
via Fm,k- At this point we note that it is impossible to 
solve Eq. 1391 analytically for general q and we will intro- 
duce an approximation justified for very long polymers. 
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Ma, Mb ^1. In this case the relevant wave vector is 
expected to be small in magnitude and consequently, the 
"end factors" as well as the wave dependence of (j)k,k' are 
negligible. Therefore, we replace them by their zeroth 
order approximations in q, 



$k.k' = -Ikh'idk + dk')\uJk X (40) 



and 



exp (cr„,feiiq • l„iLdrn) = 1 

cos (izq • /fcd)fc) = 1 (41) 

Then the bifurcation equation becomes 



,i(q,aj„i) = / pk' ,i{q_,i2i')(i>m,k' {u}7n,oj')' 

^ k' 

k^m i-' •' 



(4' 



(42) 



where again as in the case of the nematic solutions, 
4>k,k' {ijJ, ijj') has the convenient property that it maps Pj 
on Pj . Then the only mode for which the second term on 
the right hand side of Eq. 021 survives (and we have wave 
dependence) is for Pg. (For j ^ we simply re-obtain 
the nematic bifurcation equation, Ea. lHOh Consequently, 
integrating both sides over , we obtain 

c^H?) = E ^".^-(9) E ^^^'^'('^^ + dk'^oc^i'Hq) 

k k' 

(43) 

where we have defined Cm\q) ~ J dLL!mPm,ii'iT'^m) and 
where 



duj(j>k,k'{<^,'^') = -ikh'idk +dk')so, (44) 



with sq = tt"^. The rest of the analysis is similar to the 
nematic case: again we have an M x M eigenvalue equa- 
tion and we make use of the property of the geomet- 
ric factor that it depends on the types involved and not 
on the segment labels, hence J^k' — Sr' X^/c'gt' '^'^^ 
r' = A,B. Defining c^°\q) = {\/Mr)Y.ker£k\(l) and 



F^y = (l/M,M,OE™erEfc'er'^™,fc: Eq.Hbecomes 



T 



— E ^-.-'(«) E Mr'Mr"lr'lr"{dr' + d.")4'' (?) 

(45) 

Rewriting in terms of dimensionless quantities, we obtain 



co('7) 



TTUM 



+ M){1 + MPd) 



.^Fiq)GoCo{q) (46) 



with 
Go 



iM;(l + d) 



^Ml{l + d) h'PPd 



JO) 



and Co(q) = ( Yo) I (47) 



The elements of F{q) are 

fA,A^^(MA + 



1 - jo(g'A) 



{Ma - 1) 



1 - joiqlA) 



Fa,b = Fb.a - 



1 l-UoiqlA)r'^ 



Ma 



X 



1 - Jo(9'a) 

1 l-{jo{qlB)V'- 

Mb l-joiqh) 



(48) 



and 



Fb,b^^(Mb+ 



1 - Jo('7^b) 



(Mb 



JojqlB) - Uo{qlB)V'^ 
1 -joiqlB) 



Again there are two solutions for this 2x2 eigenvalue 
problem but this time the plus sign (see again Eq. I34|) 
yields the physical bifurcation density, rimps, for the mi- 
croseparated phase (mps), 

4(1 + M)(1 + mPJ) 
""""P^ " ^M 



tr(F((?)Go) + A/tr^(F(g)Go) -4dctF(g) dot G( 



/(2detF(q)detGo). (49) 

Apart from the approximations made in formulating the 
model, Sec. Eqs. 00] and constitute the only two 
further approximations. From Eq. 0^ it is observed di- 
rectly that the spinodal density of the microseparated 
phase scales with 1/M, contrary to the nematic spin- 
odal, Eq. which does not depend on M in this repre- 
sentation. Consequently, for long enough polymers the 
system will always become unstable with respect to the 
microseparated phase. Furthermore, we note that for 
infinitely long chains (M oo) the approximations be- 
come exact (and the density needs to be rescaled, hM). 
If the chains are not long, the approximations, Eqs. I4UI 
and will not be valid. An interesting case are e.g. 
rod-coil copolymers where Ma = 1 and Mb is large. The 
type-A rods will tend to form a smectic which the type-B 
tails are likely to stabilize [sol Isij. In this case, Eg. 1391 
has to be solved numerically or in some other (approxi- 
mate) way. Moreover, the ordering of the type-A rods is 
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then likely to be dominated by an orientationally ordered 
density fluctuation, e.g. possibly exp[iq • r]P2{q • cD), in- 
stead of the simple exp[iq-r] which we have in the present 
case. Finally, we note that the specification of the geom- 
etry is contained in the matrix F(q). Using other geome- 
tries, e.g. ABABAB... repeating multiblock copolymers 
or branched geometries, do not change Eos. 1391 orl49lbut 
only the form of F{q) (the only requirement is that there 
are no closed loops within the polymers ^,32i]). 



V. THE GAUSSIAN LIMIT 

In this section we will construct a consistent limit for 
infinitely long chains of our model. There are several rea- 
sons for this approach. First of all, there is a large body 
of literature dealing with so-called Gaussian chains, i.e. 
polymers which are coarse-grained on the level of the 
radius of gyration, and we want to make contact with 
those treatments P, |33| ■ Secondly, we do not fully con- 
trol the quality of the approximations, Eqs. ^U] and 1411 
made for chains of finite length. It is clear, however, 
that these approximations become exact for infinitely 
long polymers. Finally, by introducing this limiting case 
the number of effective model parameters is reduced, re- 
sulting in a conceptually simpler system. The limit of 
Ma, Mb — > oo does require that some of the other param- 
eters be rescaled as well. Additionally, we want to take 
this limit in such a way that the nematic and microsep- 
arated bifurcation densities remain of the same order of 
magnitude so that we can compare them. This extra 
requirement is non-trivial as can be seen from Eqs. 1351 
and 1491 because "hmps scales with 1/M and thus vanishes 
for long polymers. We can cure this divergence in a 
somewhat unconventional way by letting the difference in 
thickness of the A and B segments vanish, d ^ 1. In this 
way, the incentive for MPS is much reduced and rimps 
"pulled up" to nonzero densities comparable to nnom- 
Summarizing, we take the limits 



defined as 



Ma 



— > and d ^ 1 



(50) 



whilst MpJ.\ and Ma(1 — d)^ remain finite. Furthermore, 
we keep the ratios M and I fixed, such that the type-B 
segments are subject to the same limit. Next, in order 
for the Onsager approximation to still be valid, ^a needs 
to remain much smaller than Za and therefore needs to 
go to zero even faster. This is corrected by letting the 
number density of chains go to infinity in order to keep 
total strength of the interaction, i.e. the total excluded 
volume constant. So additionally we have 



cZa — and n — > oo 



(51) 



with 2nMfiJ.\dp^ and therefore also n = 2n{MfiJ\dA + 
MB^IdB) finite. 

In the Gaussian limit, the relevant length scale is the 
radius of gyration or equivalently, the mean-square end- 
to-end distance. The mean-square end-to-end distance is 



fc.fe' 



(52) 



where <> denotes the average over a single chain |55l |. 
In a freely-jointed chain there is no orientational cor- 
relation between the segments so for our block copoly- 
mers, the mean-square end-to-end distance is simply 
= MpJ,\ + M^l^. This allows us to define the di- 
mensionless wavenumber as q — qx. 

Our reduced model has three parameters, M , I govern- 
ing the composition and the relative size of the copoly- 
meric blocks and A = Ma(1 — d)'^ describing the remain- 
ing difference in thickness between the two components 
and hence effectively setting the incentive for demixing. 

In the Gaussian limit, the determinant of G2 goes to 
zero, detG2 = —jMP{l — d)^ 0. Consequently, we 
can expand Ea l35l for small det G2 and we obtain for the 
nematic bifurcation density in the Gaussian limit. 



32(1 + mP 



32 



7rtrG2(o? =1) 71- 



(53) 



which, conveniently, is a constant independent on the 
model parameters. Setting the first element of the eigen- 
vector to one, Cnem — (l,c„cm), this is very simple in 
the Gaussian limit, Cnom = I- Therefore, at the nematic 
bifurcation the B segments are I times more strongly ori- 
entationally ordered than the A segments. 

Concerning MPS, we first calculate the elements of F 
in the Gaussian limit, 



Fa,a = 



^|i-4(> 



FA,B==FB.A = 4(l-e-^iMA(i_ 
9a ^ ^ 9b ^ 



6_ 

9I 



(54) 



Fb,b — 



96 I 9b ^ 



with 



9a = 



1 + MP 



and q^ = 



1 + MP 



(55) 



The determinant of Gq also goes to zero, det Gq = 
-j]VPP{l - df 0. Next, expanding Eq.EHIfor smaU 
det Go as well, we obtain for the bifurcation density of 
MPS in the Gaussian limit. 



^mps 



— lim 



+ M){1 + MPd) tr(FGo) 



ttM 
16(1 + M)(1 



det F det Go 
mP) tr(FGo) 



7rA2M2Z2 



detF 



(56) 
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FIG. 2: Bifurcation density for the microseparated phase vs. 
the magnitude of the wave vector for I — 1 and A = 4 and 
M = {5, 4, 3, 2, 1} (from top to bottom). The nematic bifur- 
cation density fincm = 32/7r « 10 and has no wave dependence 
but is drawn as a straight line for comparison.. Due to sym- 
metry ({M,i} {l/A/,1//}) the curves of M are the same 
for 1/M. 



with Go = limj^-i^ Gq, 



Go 



1 Ml 

Ml IVPP 



(57) 



Additionally, we note the symmetry in the A and B 
types, i.e. the following transformation {A,M,l} 
{A, 1/M, 1/1} leaves the results unchanged. Again, writ- 
ing the eigenvector as follows Cmps = (1, Cmps) we obtain 
a simple expression in the Gaussian limit Cmps = —1/ Ml. 
This is the relative order of component B over A at 
bifurcation. The minus sign is due to the difference 
of TT in phase between the density waves of A and B, 
i.e. where the density of A is enhanced the density of 
B is depressed (e''^ — —1). The absolute value 1/Ml 
is ratio of amplitudes of the two waves. The matrix 
F contains the correlations within the polymer and is 
seen to feature the so-called Debije functions, goix) = 
{2/x){l — (l/a;)(l — exp[— x])) reflecting the Gaussian 
character of the correlations. In the Leibler approach [s^ 
these appear in a similar way and therefore, the correla- 
tions are treated on the same level. 



VI. RESULTS 

A. Bifurcation Density 

In Fig. 12 we have plotted the analytical bifurcation 
density of the microseparated phase, Eq. H56|l as a func- 
tion of the wave vector for various values of M. Most im- 
portantly, all curves have a minimum for a certain wave 
vector. Interpreting the bifurcation point as the spinodal, 
where the isotropic fluid phase changes from being stable 
to unstable, the system becomes first unstable for fluctu- 
ations with a wave length corresponding to the minimum 
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FIG. 3: The minimum bifurcation density for the microsepa- 
rated phase vs. logM for A = 4 and [ = {0.25,0.5, 1,2,4} 
(right to left). The nematic bifurcation density nnem = 
32/7r « 10 is constant and drawn as a straight line for compar- 
ison. Inset: the wave length for which the bifurcation density 
of the microseparated phase is a minimum, Amin = 27r/gmin vs. 
M for the same parameters; A = 4 and / = {0.25, 0.5, 1, 2, 4} 
(right to left). 



density. We have also plotted the nematic bifurcation 
density, being a constant independent of the wavenum- 
ber q, in Fig.|2] For the curves which lie totally above the 
horizontal line, the system becomes unstable with respect 
to the nematic phase at the density fi — rinom — 32/7r. 
For a curve of which the minimum reaches below the hor- 
izontal line, the system becomes unstable with respect to 
a microseparated phase with wave length Amin — 27r/g,nin 
at the minimum density n = fi^^s^ . In Fig. |21 we have 
set the A and B segments to equal length, I — 1 and 
the demixing parameter is A = 4 . Starting with an 
asymmetric polymer, M = 5, MPS only occurs for high 
densities. Making the polymer more symmetric and de- 
creasing M to one, the curves shift to lower densities until 
at M = 1 it is at its lowest position. Upon a further de- 
crease M following the sequence M = {l,i,i,i,i}, we 
again follow the same curves in Fig.[21due to the symme- 
try (A, Af, 1} (A, 1/M, l/l] and the choice [ = 1, but 
now from the bottom to the top. 

We have numerically determined the minimum of the 
MPS bifurcation density with respect to the wave vector, 
Eq. H56|) and plotted that in Fig. O as a function of M 
for a few different I. We observe the same trend we saw 
in Fig. 12 for very asymmetric polymers, M <C 1, the 
minimum MPS bifurcation density is very high. Increas- 
ing M, the bifurcation density goes down until a certain 
value M (depending on /) after which it goes up again. 
As shown in Fig. |21 some of the curves reach below the 
horizontal line marking the stability limit of the isotropic 
phase towards nematic ordering. Consequently, in the in- 
termediate region the microseparated phase is probably 
the most stable phase, whereas for the more asymmetric 
polymers MPS is likely to be preempted by the nematic 
phase. Furthermore, there is also a dependence on Z, i.e. 
increasing the asymmetry between the A and B segments. 
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the curves shift to higher densities. Again, we note that 
the two curves for I — 0.5 and I = 2 can be mapped onto 
each other due to symmetry in the model parameters. In 
the inset of Fig. O we have plotted the value of the wave 
length A = 27r/g corresponding to h'i^s^ vs. M. There 
is a rough correspondence as a function of M in that the 
lower the MPS bifurcation densities in Fig. |31 connect to 
the higher wave lengths in Fig. I^Jinset). In general, we 
have observed that the wave length for which the MPS is 
the stable phase (over the nematic) roughly lie between 
1 and 1.5 times the mean end-to-end distance x, i.e. the 
polymers get somewhat stretched at the phase transition. 



"T — ' — I — ' — r 




B. Phase Diagrams 

Figs.01and[Slpresent the phase diagrams. We have nu- 
merically computed the model parameters for which the 
minimum MPS bifurcation density equals the nematic bi- 
furcation density. In Fig. 0] the phase diagram is given 
in terms of M vs. A for equal length segments, 1 = 1. 
For low A the incentive for MPS is too weak and the 
MPS bifurcation densities are higher than the nematic 
ones everywhere. Increasing A, the MPS becomes stable 
for M = 1 (totally symmetric diblock copolymer) and in- 
creasing A further the range of M for stable MPS grows 
correspondingly. This is not surprising as the MPS bifur- 
cation density scales simply with 1/A. The inset of Fig.0| 
shows the vertical scale logarithmically to show the sym- 
metry with respect to M 1/M. In Fig. O the phase 
diagram is plotted for M vs. I. The same observation 
as in Sec. IVI Al can be made: for asymmetric polymers, 
the nematic phase is the most stable whereas for more 
symmetric ones the MPS can be stable. Of course the 
amount of area in Fig. [S] depends sensitively on A. Note 
that I plays a very similar role as M. Naively, one might 
expect that a difference in lengths of the segments would 
also increase the tendency to microphase separate or at 
least not counteract to it. However, this is not the case, 
and only the difference in thickness, even though only 
infinitesimally small in the Gaussian limit, drives the oc- 
curence of MPS, in line with earlier work on binary mix- 
tures of rods 1221 . Potentially, length differences between 
the component rods could drive MPS within the nematic 
phase, but probing this would require the numerical so- 
lutions to the full self-consistency problem, currently be- 
yond our scope. 



C. The Density Shift along the Polymer 

The elements of the eigenvectors at the bifurcation 
as discussed in Sees. IIVBI and IIVCI contain informa- 
tion about the relative amplitude of the nascent order- 
ing with respect to the the homogeneous and isotropic 
parent phase. However, by construction these quantities 
were averaged over all segments either of type A or B. In 



FIG. 4: Phase diagram, M vs. A for / = 1. For the re- 
gion marked with "Nem", the lowest bifurcation density is 
the nematic and for the region marked with "MPS" this is 
the microseparated phase. The inset is the same phase dia- 
gram except that the vertical axis is logarithrnic to show the 
symmetry with respect to {M,l} —f {1/AI, l/l}. 
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FIG. 5: Phase diagram, M vs. I for A = 4. For the region 
marked with "Nem" , the lowest bifurcation density is the ne- 
matic and for the region marked with "MPS" this is the mi- 
croseparated phase. The inset is the same phase diagram 
except that the axes are logarithmic to show the symmetry 
with respect to {M, 1} {1/Af, l/l}. 



case of the nematic ordering, this also coincides exactly 
with the order of each of the segments individually as 
there is no orientational coupling between the segments 
and these therefore behave as being independent. How- 
ever, in case of MPS, there clearly is a spatial coupling 
between the segments and, consequently, one would ex- 
pect a different degree of ordering e.g. for segments which 
are close to the free end and those which are close to the 
joint. Those close the joint are be subjected to two coun- 
teracting density waves and will order less than those at 
the free ends. In order to to quantify these effects we 
have to compute the components of the M-dimensional 
vector Cm' (Eq. (021 )• In appendix O we explain how 
these are obtained from the type-averaged 2-dimensional 
eigenvectors by means of an additional quantity: the half 
type- averaged matrix F'. In the Gaussian limit, this M- 
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FIG. 6: Relative order along the polymer at bifurcation in 
the microseparated phase, c'}^'' (s) for s £ [0, 1/(1 + M)] and 
c'^°\s) for s e [1/(1 + M), 1]. Parameters are A = 4, [ = 1 
and M = {1,1.5,2,3,4,5,7,10} (increasing in the direction 
of the arrow). The normalization is such that the averages 
over (s G A) and Cg"' (s £ B) equal c^' — 1 and Cg^' = 
~1/{MI) respectively. The full circles indicate the "joints" of 
the A and B parts at s = 1/(1 + M). 

dimensional vector reduces to the following 2-dimensional 
eigenvector (with a prime), 

which now depends, on the continuous label s G [0,1], 
where s e [0, j:^] implies s e A and s G [ ^^^^ , 1] im- 
plies s € B. In Figs. El and 13 we plot the components of 
Cq(s) along the polymer (as a function of s) for increasing 
M and I respectively. The demixing parameter is taken 
to be A = 4. 

In Fig. El we start from the symmetric case, M — 1 and 
I = 1 where the profile is also symmetric around s = 0.5. 
All A segments have positive order and all B segments 
have negative order and the average of A and B is +1 and 
— 1 respectively as expected. Increasing M, the B part 
of the polymer becomes larger than the A part and the 
joint shifts to the left. The normalization remains such 
that average order of the A segments is still 1 and that 
of the B segments is —1/M. However, it is remarkable 
that the B segments close to the joint obtain a positive 
order with increasing M , i.e. they order with respect to 
the density wave of A instead of that of B. This is due 
to the fact that in the polymer there is much more ma- 
terial from the B part. Consequently, this effect becomes 
stronger for larger M. In Fig. \7\ we start again from the 
symmetric case, M = 1 and 1 = 1. Subsequently, the 
ratio of lengths I is increased and we see that the deriva- 
tive of the profile to s jumps at the joint. Furthermore, 
also here, the joint shifts to positive values and the A 
segments have a much more constant profile than the B 
segments. By increasing I while M remains constant one 
effectively increases the amount of material in the B part 
of the polymer. Therefore, it is not surprising that the 
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0.2 0.4 0.6 0.8 1 

S 

FIG. 7: Relative order along the polymer at bifurcation in 
the microseparated phase, c'j^^s) for s e [0, 1/(1 + M)] and 
c'-^°\s) for s e [1/(1 -I- M), 1]. Parameters are A = 4, M = 
1 and I = {1,1.5,2,3,4,5,7,10} (increasing in the direction 
of the arrow). The normalization is such that the averages 
over cj^\s € A) and Cg''-'(s £ B) equal c^' = 1 and c^' = 
-~1/{MI) respectively. The full circles indicate the "joints" of 
the A and B parts at s = 1/(1 -I- M). 

point of zero order shifts to the right. Additionally, the B 
segments are much longer and therefore the spatial corre- 
lations persist over larger s explaining the more smooth 
profile on the B side. It has to be noted that some of the 
profiles (especially for higher values of M and I in Figs. El 
and[7|) are taken at bifurcation densities far above the ne- 
matic bifurcation. We have nevertheless included them, 
being instructive in explaining the observed trends. 

VII. CONCLUSION 

We have considered a fluid of freely-jointed hard di- 
block copolymers. The two polymer blocks A and B 
consist of slender Onsager rods of different dimensions 
interacting via hard body repulsion only. We apply a 
DFT approach in the second virial approximation from 
first principles, and analytically construct local solutions 
to the stationarity equations, by means of a stability (bi- 
furcation) analysis of the isotropic phase. Spatial as well 
as orientational degrees of freedom are taken into account 
and consequently we obtain the spinodal densities for 
both the microseparated and the nematic phases. It is 
shown that for long polymers the system always becomes 
unstable with respect to the microseparated phase first. 
Consequently, this means that entropy can induce MPS 
in much the same way as it has been found to induce 
other forms of spontaneous ordering before. Further- 
more, the mechanism is determined solely by the (dif- 
ference in) dimensions of the rods and therefore has a 
conceptually simple geometric origin. 

In order to make contact with the literature on ther- 
motropic block copolymers we take the limit of infinitely 
long polymers in which the approximations become ex- 
act. In addition, by assuming a vanishing difference in 
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thickness of the two types of rods, we can still study 
the competition of the microseparated with the nematic 
phase. We present phase diagrams in terms of model pa- 
rameters showing the regions of stable microseparated or 
nematic ordering. We also present the order along the 
polymer at the bifurcation of the microseparated phase. 

In the present study, we have solved the stationarity 
equations up to first-order in a bifurcation analysis. This 
yields, apart from the location of the spinodal or bifur- 
cation density, only the magnitude of the density wave 
vector and the spherical harmonic mode to which the 
isotropic solution becomes unstable. However, the sym- 
metry of the bifurcating microseparated solution is typ- 
ically determined by one or more mutually independent 
(but equally long) vectors spanning the periodic phase 
(e.g. lamellar, hexagonal or bcc). In order to obtain in- 
formation on the mutual orientation of these lattice vec- 
tors, and thus on the symmetry of the phase, a higher or- 
der bifurcation analysis should be performed 31, 56, 57] . 
From these higher order bifurcation equations, it is also 
possible to determine whether the phase transition is of 
first or second order and in the latter case one could in 
principle go on to construct the full equilibrium solution 
far away from the bifurcation point j3l| . 

We have not checked the validity of the approxima- 
tions, Eas. I4UI and 1411 for finite values of M. However, 
we can make a crude estimate, a posteriori, by con- 
cluding from Fig. that the bifurcating wave length is 
of the order of the mean square end-to-end distance, 
A = A/x ~ 1. Consequently, the wave vector is ap- 
proximately, q = 27r/A 27r and if we assume for a 
moment that the type-A rods and type-B rods are more 
or less equally long, then the mean-square end-to-end 
distance is ~ Ml\. This in turn implies that the 
next order corrections in Eg s . 1401 and 14 1 1 will be of order 
{\qlAf ~ (i X 2-nflM - 10/A/. (In fact the first order 
correction in (^^^a) in Eq.^Jdoes not contribute to the 
value of the bifurcation density, but only to the form of 
the eigenfunction.) Consequently, already for this crude 
test case, the length of the polymer should be at least 
longer than 10 {M > 10) in order for the corrections 
to be smaller than the leading term. This suggests that 
much higher values of M are required for the present ap- 
proach to yield quantitative agreement with the "true" 
behaviour. 

In any case, it would be very interesting to extend the 
present approach to finite values of AI. However, this is 
not straightforward, as the correlations within the chain 
would become non-Gaussian. One strategy could be to 
solve Eq. directly numerically but this could become 



tedious for large numbers of segments. Another strat- 
egy would be to make an expansion in 1/M using the 
Gaussian limit as a reference state. This last route was 
followed by Fredrickson and Hclfand for Leibler's di- 
block copolymers and the results were confirmed by sim- 
ulations [33 • Indeed, there is a need for such a better- 
than-Gaussian treatment, especially when the typical or- 
dering length scales are of the same sizes as the com- 
ponents, e.g. for side chain liquid crystalline polymers 
forming a smectic |5^ [S^ . 

As already mentioned in the introduction, there is 
as yet no experimental system exhibiting MPS due 
to the mechanism described in this paper. However, 
considering the ongoing progress in the field of bio- 
engineering [2^ , it may become possible to prepare 
such a system. We mention again the possibility of long 
and thin polymers connected to TMV rods in an appro- 
priate solvent. The solvent may be a problem as we have 
the double requirement that the polymers are at their 9- 
point and that at the same time the TMV rods still act 
as hard particles. Still, such a system of entropic rod-coil 
copolymers could be directly compared to the simulation 
studies of Refs. 50, 51]. Additionally, it would be de- 
scribed by Eq. 12^1 which would than have to be solved 
for the case of Ma — 1 and Mb large. In a more gen- 
eral context, it becomes increasingly clear that entropy- 
induced effects play a prominent role in vivo [bH ]. and it 
may be that similar mechanisms as described here pre- 
vent demixing tendencies due to local constraints [60j. 
On the other hand, the mechanism may also be of rele- 
vance in thermotropic systems where the two components 
of block copolymers also have short-range anisotropic re- 
pulsions which are usually of different range. In any 
case, observing entropy-induced microphase separation 
in monodisperse systems would certainly be an interest- 
ing experimental challenge. 
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APPENDIX A: THE FOURIER TRANSFORMED 
SEGMENT-SEGMENT MAYER FUNCTION 

The Mayer function (j)k,k' of two cylindrical rodlike seg- 
ments k (with dimensions Ik, dk and coordinates (r/j, uJk)) 
and k' (with ly , dk' and (r^/ , cJk')) interacting via a hard- 
core potential (i.e. = oo, if overlap/no overlap) is given 
by 



I if overlap 
if no overlap 



(Al) 



We decompose the spatial vector rk,k' 
of the orientations, 



Tfc — Ffc/ in terms 



rfc,fe' = XkUJk + Xk'iOk' + Xk,k'l^k,k' 



(A2) 



with u)k.k' = {i^k X '^k')/\i^k X i^k'l the unit vector in the 
perpendicular direction. There is overlap between the 
two rods for the following ranges of the coefficients, Xk S 

[-/fe/2,/fe/2], Xk' € [-lk'/2,lk'/2] and Xk^k' e [-{dk + 
dk')/2, {dk -\- dk')/2]. Next, the Fourier transform of the 
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Mayer function (f>k.k' is given by 



(A3) 

where the volume of the infinitesimal element is given by 
drk,k' = \i^k X LOk'\dxkdxk'dxk,k> ■ Consequently, 



with the asterisk denoting the complex conjugate and z 
some unit vector. It is now directly seen that the Legen- 
dre polynomials are eigenfunctions of (j)k,k'{'^ • 



du}'^k,k'{i^ ■ i^')P]{'^' ■ z) = -lkh'{dk +dk')sjPj{Lb ■ z). 



(B5) 



(/)fc,fe'(q,d;fc,a)fc') = 














)/2 


- \ilik y.ujk'\ / dxk 


/ dxk' 




dxk.k' X 








0/2 



exp[-i(a;fcq • uJk + Xk'q ■ ujk' + a^fc^fc-q • ujk.k')] 
= -Iklk'idk + dk')\uJk X uJk>\jo (I'fcq ■ i^k) X 

jo {^h'fl ■ Ld'k,) jo {^{dk + 4')q • <^k,k') (A4) 

with the spherical Bessel function of zeroth order given 
by jo{x) — sinx/x. In the Onsager limit of very slen- 
der rods, IkJk' > dk,dk' while hh'idk + dk') stays fi- 
nite. In our system, we expect the wave length of the 
microseparated phase to be at least of the order of the 
lengths of the segments (although for large number of 
segments it is even much larger). Consequently, in this 
case, \(dk + dk>)q\ ^ 1 and we use the leading order, 
which is jo{^{dk + dfeOl ' '^k,k') — 1- Then, our final 
result for the Mayer function is 



(^fc,fc'(q,d;fe,a)fe') = -hh'idk + dfcOI'^fe x Lo'k'\'^ 
jo i^hq. ■ <^k) jo i^h'q. ■ uj'k') 



APPENDIX B: THE EIGENFUNCTIONS OF 

4,fc/(ii-cl') FOR q = 

For q = 0, the Fourier transformed Mayer function is 

4>k,k'{'^ ■ i^') = -hh'idk + dk')\Lo X Lu'\ (Bl) 
= -hh' idk + dk' ) Vl-iio-oj')^ (B2) 

and is therefore uniaxial, i.e. dependent on a single planar 
angle 7 = arccos(a) • a)'). Therefore, we can expand it in 
terms of Legendre polynomials 



ik,k' (ti • a)') = -hh' idk + 4' ) X! 

3=0 



2j + l 
47r 



S-jPjiLU ■ LU'), 

(B3) 

with Sj — 2tt dx\/l — x-^Pjix). Then, using the de- 
composition in terms of spherical harmonics Yj^i, we can 
rewrite this as 



4>k,k' (w • to') = -hh' idk + dk')x 
00 j 

EE 

j=0 %=-3 



'^-l±ls,Y,,{u;-z)Yl,iCo' -z), (B4) 



APPENDIX C: SPATIAL ORDER WITHIN THE 
POLYMER 

It is possible to calculate the bifurcating order within 

the polymer. In case of freely-jointed chains in the ne- 

(2) 

matic phase this is trivial as this exactly Cr for a seg- 
ment of type T. However, in case of MPS, segments of 
type A close to the "joint" with B segments will typi- 
cally be more affected by the B part of the polymer than 
segments of type A far away from the joint. This order 
within the polymer can be obtained by calculating the 
elements of the M-dimensional vector Cq with elements 

c^m and TO e {I,-- - ,M} (see Eq. ESJ- Therefore we 
proceed by defining the matrix F' (with a prime) 

Pm,k', (CI) 



F' 



1 



E 

fc'er' 



where the average is only performed over the second la- 
bel and therefore ^, is M x 2 dimensional. Then, if 
the bifurcation density for the microseparated phase rimps 

(A5) 

and the corresponding eigenvector Cmps has been calcu- 
lated beforehand (from Eg. 149(1 . Cq can be computed by 
evaluating 



- J' 'jqC 



Ail + M){1 + MP d) 
The elements of F' are given by 



mps ■ 



(C2) 



F' A A 




2-(j0(g/A))'"'^-(j0(g/A))^^"'" ^ 

1 - jo('7'A) 



p, _ 1 , , ,.Mj,-m 1 - (jo(g^B)) 

^™eA3-]^Oo(«^A)) ^ 



J\/b 



(C3) 



(C4) 



p, _ 1 {nl ^V»-Ma-1 1 - (jo(g^A)) 

Ma l-jo(g«A) 



Ma 



(C5) 



, 2-(jo(<?/b)) 

Mb 



m-l-AfA 



1 



ijoillB)) 



1 -io(g^B) 



(C6) 
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where m € {1, • • • , Ma} when m G A and m € {Ma + 
1, • • • , M} when m G B. For each of these elements again 
holds that the average of m yields the matrix F (see above 
Eq. 133 and Eq.gHIl, i.e. 



(C7) 



^B,B(seB) = - 2 
9b V 



exp 



exp 



9i 



(1 



l + Af 
M 

. + M 



M 



(Cll) 



In the Gaussian limit, we have to define a continuous 
"label", s = m/M, with m and M going to infinity such 
that s keeps its value. Consequently, s G [0, 1] and F' 
becomes 



FkAs e A) 



6 , „ 

-2- 1 2 - exp 



9i 

6 



s(l + M) 



exp 



- 1^(1 -.(l + Af)) 



(C8) 



where s £ [0, j^^] when s & A and s € [jq^j 1] when 
s S B. Note that in the Gaussian limit F' is simply a 
2x2 matrix, however, with s-dependence. Consequently, 
unlike F, F' is not symmetric. And additionally, also the 
M-dimensional eigenvector becomes 2-dimensional, 



n,B(se A) 



9i 



1 — exp 



6 



exp 



-^(1 
6 ^ 



s(l + M)) 



(C9) 



AO) 



cr'is £ A) 



(C12) 



F^Ai^ e B) 



9i 



1 — exp 



6 



exp 



6 ^ M 



M 



M' 



Finally, it has to be noted that in the Gaussian limit, 
(CIO) first the product of Gq and c,„ps has to be taken and 
only then the limit can be applied to (GoCmps). 



